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BY 


E  A  SKELTON 


Summary 


'  The  Green's  functions  in  the  region  between  planes  on  which 
either  Neumann  or  Dirichlet  boundary  conditions  apply  are  obtain¬ 
ed  in  the  form  of  double  Fourier  series  expansions.  A  closed 
form  expression  for  the  sound  field  in  an  axisymmetric  water 
filled  channel  is  obtained  in  order  to  illustrate  the  usefulness 
of  the  expansions.  Problems  associated  with  numerical  evaluation 
are  discussed  and  some  numerical  examples  are  given. 
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1.  INTRODUCTION 


TZTrrr- 


The  expansions  of  the  free-space  Green's  functions  of  the 
reduced  wave  equation  in  Cartesian,  cylindrical  and  spherical 
polar  coordinate  systems  have  been  used  extensively  in  the 
literature  to  solve  boundary  value  problems,  in  fluid-structure 
interaction,  in  which  at  least  one  of  the  coordinate  directions 
is  infinite.  One  of  the  Green's  function  expansions  tabulated  by 
Skelton  CID  appears  to  be  relatively  new  to  the  literature  (see 
equation  (2.2)  of  this  memorandum).  James  C2]  and  Skelton  C3J 
have  used  it  to  facilitate  solution  of  the  problem  of  infinite 
pipes  and  layered  cylinders,  respectively,  that  are  excited  by  an 
interior  point  source  of  sound:  Fuller  has  used  it  to  investigate 
infinite  shell  vibration  C43  and  the  sound  transmission  into  an 
idealised  model  of  an  aircraft  fuselage  C5D.  Thus,  when  new  forms 
of  Green's  functions  become  available  they  are  invariably 
accompanied  by  the  solution  of  a  class  of  problems. 

It  is  desirable  to  extend  the  range  of  Green's  function 
expansions  to  enable  solution  of  problems  involving  finite 
geometries.  For  example,  the  sound  field  inside  a  finite 
cylindrical  enclosure,  and  the  radiation  of  sound  from  the 
elastic  vibrations  of  its  surface  are  problems  of  considerable 
interest  in  the  field  of  noise  and  vibration.  Solution  of  these 
problems  would  be  facilitated  if  the  Green's  function  in  the 
region  between  planes,  on  which  Neumann  or  Dirichlet  boundary 
conditions  are  applied,  were  available  as  double  Fourier  series 
expansions  in  the  axial  and  circumferential  coordinates  of  a 
cylindrical  coordinate  system.  The  Green's  function  expansions 
when  the  source  is  on  the  z-axis  are  well  known  from  sound 
propagation  work  C63,  but  the  extensions  to  an  off-axis  source 
have  not  been  found  in  the  literature. 


2.  GREEN'S  FUNCTIONS 


( a )  General 

An  infinite  slab  of  acoustic  fluid  is  located  between  the 
planes  z=0  and  z=h  cn  which  either  Neumann  (hard)  or  Dirichlet 
(pressure  release)  boundary  conditions  apply.  A  time-harmonic 
point  source  of  sound  of  free-field  pressure  exp( ik| R-R^ I ) / I R-R^ I 

is  located  in  the  acoustic  fluid  at  the  cylindrical  coordinates 
(r,(p,2)  =  (r^,0,2^).  k=u)/c  is  the  acoustic  wavenumber  in  which 

u)  is  the  angular  frequency  and  c  is  the  sound  velocity.  The 
geometry  is  shown  in  Figure  1.  The  time-harmonic  factor, 
exp(-lu)t),  is  omitted  from  all  equations. 

The  pressure  in  the  fluid  layer  satisfies  the  inhomogeneous 
reduced  wave  equation 

(V^+k^ )p(R,R„)  =  -4w6(R-R  )  (2.1) 

r  _  -0  -  -0 

together  with  the  selected  Neumann  or  Dirichlet  boundary 
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conditions  at  z=0  and  z=h.  The  pressure  may  thus  be  written 
as  the  sum  of  the  free-space  Green's  function  Pq(R,R^)  together 

with  a  scattering  term  pg(R,R^)  which  satisfies  the  homogeneous 

form  of  equation  (2.1)  and  whose  amplitude  is  determined  by  the 
boundary  conditions. 


(b)  Free-Space  Green's  Functions 

The  free-space  Green's  function  of  the  reduced  wave  equation 
in  cylindrical  polar  coordinates  in  a  form  which  is  appropriate 
to  the  application  of  additional  boundary  conditions  on  cylind- 


rical  boundaries  is  C13 

f  " 

00 

(i/2)  y  e^cos(n((i) 

1 

J  ( Yr)H^( yr^ )exp( iaCz-z„3)da 
n  n  0  0 

n=0 

-00 

rv<r 

Pg'B'Bo’  =  ■ 

00 

O 

OQ 

(i/2)  y  e  cos(n<|)) 

f 

J  (yr  )H_(yr  )exp(  i(x,Cz-z^3)da 
non  0 

n=0 

-00 

r>,r^ 

(2.2) 

where  J_  and  H_*J_+iY_  are 
n  n  n  n 

Bessel  and  Hankel  functions. 

respectively ;  and  ,  with  Im(-Y)>/0  in  order  to  satisfy 

the  radiation  condition. 

The  infinite  series  in  equations  (2.2)  converge  only  in  a 
generalised  function  sense.  They  converge  numerically  provided 
that  dissipation  is  included  in  the  system  by  setting,  for 
example,  the  sound  velocity  c  to  be  the  complex  value  cd-in) 
where  n  is  a  very  small  hysteretic  loss  factor.  However,  when 
r=r^  the  series  cannot  be  evaluated  numerically  even  when 

dissipation  is  included. 

An  alternative  form  of  equation  (2.2)  which  is  required 
later  is  determined  from  an  application  of  a  variant  of  Graf's 
addition  theorem  for  Bessel  functions  C73,  viz.. 


00 

H  (yJCr^+r^-2rr  cos<p3)  =  )  e  H  {yr  ) J„(-Yr)cos(n4»)  r  >,r 

0  0  0  /,nn  on  o 

n=0  (2.3) 

in  which,  if  it  is  necessary  only  to  interchange  r^  and  r. 

The  required  form  is  obtained,  by  substituting  equation  (2.3) 
into  equation  (2.2),  as 

00 

=  (i/2)  [  H  (>J Cr^+r^-2rr  cos03)exp(iaCz-z  3)d(x  (2.4) 


( c )  Application  of  the  Method  of  Images 

It  is  well  icnown  C63  that  the  scattered  pressure  due  to  the 
presence  of  a  single  rigid  plane  boundary  is  equivalent  to  the 
pressure  of  a  second,  or  'inage',  point  source  which  is  located 
at  an  equal  distance  on  the  opposite  side  of  the  boundary  to  the 
original  source,  both  sources  having  equal  magnitude  and  phase. 
Similarly,  the  pressure  due  to  the  presence  of  a  single  pressure 
release  plane  boundary  is  equivalent  to  the  pressure  of  an  image 
source  of  equal  magnitude  to,  but  out  of  phase  with  the  original 
source.  For  the  case  to  be  considered  here  in  which  the  source 
is  located  between  two  parallel  plane  boundaries  the  pressure  in 
the  fluid  may  be  calculated  by  summing  an  infinite  series  of 
image  sources  obtained  as  a  result  of  successive  mirror  reflect¬ 
ions  of  the  source  in  the  boundaries.  The  locations  of  the  image 
sources  (see  Figure  1)  on  the  vertical  axis  are  straightforwardly 
obtained  as  the  four  groups  of  coordinates 

z  =  2(m+l)h+z 

0 

z  =  2<in+l)h-z^ 

0 


z  =  -2mh+z^ 
0 


z  ~  ~ 2 "Z^r 

'The  three  combinations  of  boundary  conditions  to  be 
considered  are  (i)  both  upper  and  lower  boundaries  are  rigid, 

(ii)  both  upper  and  lower  boundaries  are  pressure  release,  and 

(iii)  the  upper  boundary  is  pressure  release  and  the  lower  bound¬ 
ary  is  rigid. 


(d)  Both  Boundaries  are  Rigid 

When  both  of  the  plane  boundaries  are  rigid  it  is  clear  that 
all  the  image  sources  have  the  same  phase  as  the  real  source. 

'The  radial  and  azimuthal  coordinates  of  the  image  sources  are 
also  identical  to  those  of  the  real  source.  In  order  to  calculate 
the  Green's  function  in  this  case  it  is  sufficient  therefore  to 
replace  the  exp(-iaz^)  term  in  equations  (2.2)  and  (2.4)  by  the 

infinite  sum  I(a,z  ),  where 


I(a,z^)  =  2  Cexp( -iaCz^-2mh3)  +  exp( -iaC-Zjj-2mhD ) 
m=0 


+  exp( -iaC2 (m+1 )h+z  3)  +  exp( -iaC2(m+l )h-z  3)1  (2.5) 

0  0 


7 


This  expression  may  be  simplified,  to  give 


I(KrZ  )  -  2cos(az  )  )  exp(2i 
0  0  L 

in=-oo 


2iounh) 


(2.6) 


The  generalised  function  identity  obtained  from  the  Poisson 
summation  formula  E8D 


^  exp(  2irima/d)  "  d  ^  6(a-md) 


(2.7) 


m=-«» 


in=-«> 


which  is  valid  for  integration  purposes,  allows  equation  (2.6)  to 
be  rewritten  more  usefully  as 


I(oi,,z  )  =  (2‘ir/h)cos(az  )  )  S(a- 
0  0  L 

m=-«» 


mir/h) 


(2.8) 


When  the  expC-iotz^j)  term  in  equations  (2.2)  is  replaced  by  this 

expression  for  I(a,z^)  the  a-integration  is  straightforward 

and  the  Green's  function  of  the  reduced  wave  equation  between  two 
parallel  rigid  planes  takes  the  required  form 


00  00 

(  iit/h)  )e„cosn(|))  J^(  Y^r  )H  ( Y„r^)cos(mirz^/h)exp(  imirz/h) 
i,  n  i,  n  m  n  m  0  o  ^ 

n=0  m=-<»  , 

r,<ro 

00  00 

( iir/h)  )e„cosn4>)  J„(  Y«,r„)H„(Y_r  )cos(mirz„/h)exp(  imirz/h) 

£  n  £,  n  m  0  n  m  o  ^ 

n=0  m=-<» 

r>.ro 
(2.9) 

in  which  Yj„=(k^-m^7r^/h^  )  ^^^ .  This  may  be  simplified  to 


00 

r 


00 

r 


(i‘n/h))e  co5n<t)e  J  (>  r)H  (y  r  )cos(inir2  /h)cos(inir2/h) 
£  n  m  n  m  n  m  o  o 

n=0  m=0  - 

r^r 

o 

00  00 

(  iir/h)  >e„cosn(p)e_J  (Y_r^)H  (Y«r  )cos{mir2_/h)cos(mirz/h) 
£,  n  L  m  n  m  o  n  m  o 

n=0  ra=0  . 

r>,ro 

(2.10) 


Similarly,  replacing  the  exp(-iaz^)  term  in  equation  (2.4)  by 
the  expression  for  I(a,z  )  in  equation  (2.8),  integrating  with 
respect  to  a  and  rearranging  the  resulting  infinite  sum  gives 


an  alternative  form  of  the  Green's  function  of  the  reduced  wave 
equation  between  two  parallel  rigid  planes r  viz.. 


(  iir/h)  Ve  H  (v  J  Cr^+r^-2rr^cos^D  )cos(mirz  /h)cos{mirz/h) 
^mom  o  o  0 


m=0 


(2.11) 


( e )  Both  Boundaries  are  Pressure  Release 

When  both  of  the  plane  boundaries  are  pressure  release  it  is 
clear  that  a  phase  change  occurs  at  each  reflection.  In  order  to 
calculate  the  Green's  function  it  is  again  sufficient  to  replace 
the  exp(-iaZp)  term  in  equations  (2.2)  and  (2.4)  by  an  infinite 

sum  I(a,Zjj)  which,  in  this  case,  is 


I(a,z^)  =  2  {exp( -iaCz^-2mhD )  -  exp( -iaC-Zjj-2mh3 ) 
m=0 

+  exp( -iaC2(m+l)h+z^3)  -  exp( -iaC2(m+l )h-z^3 ) 3  (2.12) 

As  in  the  previous  case  this  may  be  simplified  considerably, 
viz.  , 


00 

I(ot,z^)  =  -2isin(az^)  ^  exp(2ioimh)  (2.13) 

m=-» 


which,  after  applying  the  generalised  function  identity,  equation 
(2.7),  may  be  written  as 

00 

I(a,Zp)  =  -(2iir/h)sin(aZjj)  ^  6(a-mir/h)  (2.14) 

m=-«» 


When  the  exp(-iaz^)  term  in  equations  (2.2)  is  replaced  by  this 
expression  for  I(a,z^)  the  a-integration  is  again  straight¬ 
forward,  and  the  Green's  function  of  the  reduced  wave  equation 
between  two  parallel  pressure  release  planes  is 


Pg^B'Bq^ 


:  ir/h)  )e„cosn(p)  J  ( Y„r  )H  ( v  r„ )  sin(m7rz,/h)exp(  imirz/h) 
L  n  n  'm  n  'm  0  o  ^  . 

n=0  m=-«» 

00  00 

ir/h))e  cosn<pyj  (y  r  )H  (y  r)sin(mirz  /h)exp(  imirz/h) 
^n  £nmonm  o  . 


n=0  m=-« 


in  which  y^=(k^-m  n^/h  )  .  This  may  be  rewritten  as 


Pq  ^  B 'Bq  ^ 


( 2iri/h) 


))e  cosn4i)J_(Y  r)H  (Y«*‘^>sin(mirz„/h)sin(mir2/h) 
in  ^i  n  'TBi  n  'm  0  o 

n=0  m=l  '  0 


00  00 

(2iri/h)Ve  cosn0)J  (Y„r  )H  (Y«r )  sin(mirz  /h)sin(mw2/h) 
in  ’^£,n  inonm  o 

n=0  in=l  '  0 

(2.16) 


Similarly,  replacing  the  exp(-iazQ)  term  in  equation  (2.4)  by 

the  expression  for  I(a,2^)  in  equation  (2.14),  integrating  with 

respect  to  ot  and  rearranging  the  resulting  infinite  sum  gives 
an  alternative  form  of  the  Green's  function  of  the  reduced  wave 
equation  between  two  parallel  pressure  release  planes,  viz.. 


Pg*-'-o'  =  (2iri/h)^H^(Vjji^Cr^+r^-2rr^cos<>3)sin(mirZjj/h)sin{niir2/h) 

m=0 


(2.17) 


( f )  One  Rigid  Boundary  and  One  Pressure  Release  Boundary 

The  case  in  which  the  boundary  at  2=0  is  rigid  and  the 
boundary  at  z=h  is  pressure  release  will  now  be  considered. 

When  one  boundary  is  rigid  and  the  other  is  pressure  release, 
phase  changes  occur  only  at  reflections  at  the  pressure  release 
boundary.  In  order  to  calculate  the  Green's  function  it  is  again 
sufficient  to  replace  the  exp(~iaz^)  term  in  equations  (2.2) 

and  (2.4)  by  an  infinite  sum  I(a,z  )  which,  in  this  case,  is 

0 


I(a,z  )  =  y  ( -1  )”'Cexp( -iaCz  -2mh3)  +  exp(-ia,C-2  -2mh3) 

0  i  ^  0  ^  0 

m=0 

-  exp( -icx,C2(m+l)h+2  3)  -  exp( -ioi.C2(m+l  )h-z  3)1  (2.18) 

o  0 

As  in  the  previous  cases  this  may  be  simplified  considerably, 
viz. , 


00 

I(a,z^)  =  2cos(az^)  ( -1  )'“exp(  2iamh)  (2.19) 

in=-oo 


This  is,  however,  a  form  which  is  not  suitable  for  a  direct 
application  of  the  generalised  function  identity,  and  it  is 
necessary  to  rewrite  equation  (2.19)  as 


00 


00 


I(a,z^)  =  2cos(oiz^)£2  ^  exp(4iounh)  -  ^  exp(2i(xinh)}  (2.! 


The  generalised  function  identity  may  now  be  applied  to  each  of 
the  infinite  sums  in  this  equation,  allowing  it  to  be  written  as 


I(a,z^)  =  (  2ir/h)  cos  (az^ )  £  ^  6  ( a-mir  /  2h )  -  ^  6(a-mir/h))  (2.21) 


which  may  itself  be  simplified  to 


I(a,z  )  =  (  2ir/h)cos(az  )  Y  S(a- (  2m+l )  ir/2h)  (2.22) 

0  0  L 


When  the  exp(-iaz^)  term  in  equations  (2.2)  is  replaced  by  this 

expression  for  I(a,z^)  the  a-integration  is  straightforward 

and  the  Green's  function  of  the  reduced  wave  equation  between  a 
pressure  release  plane  and  a  parallel  rigid  plane  is 


(  iir/h)  )e„cosn0)  J„(  )cos  (a  z^  )exp(  ia  zl 

L  n  i,  n  m  n  mo  mo  m 


n=0  m=-<» 


( iiT/h)  )e„cosn(}i)  J  (y  r^)H  ( Y_r)co5(a  z  )exp(iot  z) 
L  n  L  n  monm  mo  ^  m 


n=0  m=-<» 


(2.23) 


in  which  a  =(  2m+l  )ir/2h,  and  y  =  ( )  ^  ^  ^  .  This  may  be  re- 

m  mm 

written  as 

00  00 

( 2iir/h)  )e„cosn<J))  J  (y  r)H  (y  r  )cos((x  z  )cos(a  z) 
L  n  ^nmn  mo  mo  m 


^G^-'-o^ 


n=0  m=0 


(2iir/h))e  cosn(l>)J  (y  r  )H  (y  r)cos((x  z  )cos(a  z) 
^n  ^nmonm  mo  m 


n=0  m=0 


(2.24) 


Similarly,  replacing  the  exp(-iaz^)  term  in  equation  (2.4)  by 

the  expression  for  I(a,Zp)  in  equation  (2.22),  integrating  with 

respect  to  a  and  rearranging  the  resulting  infinite  sum  gives 
an  alternative  form  of  this  Green's  function,  viz.. 


^  w 
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p-(R,R^)  =  (2iir/h)yH„(Y«JC*’^+r^-2rr„cos4)D)cos(a  2^)cos(a  z] 
u  -  -0  i,o  i!i  oo  no  n 


(2.25) 


3.  SOUND  FIELD  IN  CHANNEL 


(a)  Mathematics 


The  use  of  the  Green's  functions  derived  in  Section  2  is 
illustrated  by  considering  the  problem  of  finding  the  pressure 
due  to  a  time-harmonic  source  which  is  located  in  an  axisymmetric 
channel  containing  water.  Figure  2  shows  the  geometry.  On  the 
free  surface  at  2=h  a  Dirichlet  boundary  condition  applies,  and 
on  the  other  boundaries  defined  by  the  surfaces  z=0,  r=a  and 
r=b,  Neumann  boundary  conditions  apply.  The  pressure  in  the 
channel  may  be  written  as  the  sum  of  the  Green's  function  of 
Section  2(f),  together  with  a  scattering  term  which  satisfies  the 
homogeneous  form  of  the  wave  equation  (2.1)  and  the  boundary 
conditions  on  the  planes  z=0  and  2=h.  Thus 


p(r,(p,2)  =  Pp(r,(|),z)  +  p-(r,(}>,z) 


(3.1) 


in  which  pQ(r,$,z)  is  given  by  equation  (2.24),  and  the 
scattering  term  is  given  by 

00  00 

Pc(r,(t),z)  =  y  e^cosn0  y  CA_^J^( Y„r )+B _ H^(Y^r)  }cos(a„z)  (3.2) 

5  L  ^  L  n  m  mn  n  m  m 

n=0  m=0 


The  unknown  coefficients  A  and  B  may  be  determined 

mn  mn 

by  imposing  the  Neumann  boundary  conditions  3p( r,0,z) /ar=0  on 
the  cylindrical  surfaces  r=a  and  r=b,  to  give  the  matrix 
relation 


fA  I 
mn 

B  ^ 
mn. 


(-2iricos(a  z„)  /hD„^) 
m  0  mn 


rJA<V„a)H  (y  r„)-| 

n  ro  n  m  n  m  n  m  o 

Lnm  nmJLnmonmJ 


where 


mn  n  m  n  ra  n  m  n  m 


(3.3) 


(3.4) 


is  the  dispersion  relation  whose  zeros  u)*tj _ are  the  natural 

mn 

frequencies  of  the  system. 


From  a  computational  point  of  view  it  is  shown  below  that 
the  alternative  form  of  the  Green's  function  given  by  equation 
(2.25)  is  better  suited  to  the  numerical  evaluation  of  equation 
(3.1) . 


(b)  Converaence  of  Series 


In  order  to  investigate  the  convergence  of  the  series  in 

equations  (2.24)  and  (3.2)  it  is  necessary  to  consider  the 

behaviour  of  the  terms  as  m-»  «.  In  this  case  y_  is  a  mono- 

in 

tonically  increasing  positive  imaginary  number  say.  The 

asymptotic  forms  of  the  Bessel  functions 


J„(ix  )  'v<  expCx  +i(nir/2+ir/4)D/(2iirx  ) 
n  m  m  m 


H  (  ix  )  '»<  2expC-x  -i  (nii/2+ir/4 )  3/ (  2iirx  ) 
n  m  m  m 

are  then  used  to  obtain  the  asymptotic  forms 


“  expCx^(r-b)+x  (r„-b)3 
mn  n  m  m  mo 


B  H  (y  r)  “  expCx  (a-r)+x  (a-r  )D 
mn  n  m  m  m  o 


(3.5) 


(3.6) 


(3.7) 


of  the  coefficients  in  the  series  expansions  of  the  scattered 

pressure,  p™ ( r ,0, 2 ) .  Clearly,  unless  r=r  =b,  equation  (3.6) 
o  0 

decreases  exponentially  as  Xj^-»  <».  Additionally,  unless  r=rp=a, 

equation  (3.7)  decreases  exponentially  as  ».  Thus,  unless 

both  source  and  receiver  are  situated  on  the  inner  or  the  outer 
cylindrical  boundary,  rapid  convergence  of  the  scattered  pressure 
summation  is  assured. 

If  the  same  reasoning  is  applied  to  the  terms  in  the 
Green's  function  expansion  given  by  equation  (2.24)  it  becomes 
evident  that 


)  «  expCx^ ( r " r  )  3 , 
nm  n  mo  m  o 


*  eXpCX„  (  r„ -r  )  3  , 
nmonm  mo 


(3.B) 


(3.9) 


Thus  unless  convergence  is  assured  but  at  a  rate  which  is 

less  than  that  of  the  scattered  term.  When  r=r  the  series 

0 

does  not  converge  numerically  -  it  exists  only  as  a  generalised 
function.  The  terms  in  the  alternative  form  of  the  Green's 
function  decay  exponentially  as  «»  unless  both  and 

♦=0.  Thus  the  alternative  form  of  the  Green's  function  is 
suitable  for  numerical  evaluation  everywhere  except  when  the 
receiver  is  on  a  vertical  line  through  the  source. 

The  results  of  this  Section  are  applicable  to  other  problems 
of  this  kind.  The  main  point  that  arises  is  that  the  alternative 
forms  of  the  Green's  functions  should  be  used,  where  practical, 
when  numerical  computations  are  required. 
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4.  NUMERICAL  RESULTS 


A  Fortran  program  has  been  written  to  calculate  the  sound 
field  due  to  a  time  harmonic  source  located  in  an  axis3niinietric 
channel.  The  first  version  of  the  program  computed  the  Green's 
function  using  equation  (2.24).  Severe  numerical  problems  were 
encountered.  When  the  alternative  form  of  the  Green's  function 
equation  (2.25)  was  used  in  the  computations  rapid  convergence 
was  achieved,  thus  helping  to  validate  the  analysis  of  the 
previous  Section. 

In  the  plots  of  Figure  3  the  following  constants  in  SI  units 
were  used: 


h=0.5 

a=1.0 

ro  =  l .  5 

z  =0.25 
0 


c=1500.0 
b=2.0 
r=l .  5 
2=0.25 


Thus  both  source  and  receiver  are  located  at  midwater.  The  sound 
levels  are  in  dB  ref.  1  micropascal,  and  are  to  be  referred  to  a 
free  field  source  level  of  120dB  at  Im. 

In  Figure  3A  the  angular  separation  between  source  and 
receiver  is  (ti  =  90°.  The  sound  level  spectrum  consists  of 
numerical  values  at  8Hz  spacing,  which  are  joined  by  straight 
lines.  Below  the  cut-on  frequency  of  the  first  mode  at  750Hz  the 
sound  pressure  is  negligible.  Above  this  cut-on  frequency  the 
spectrum  is  dominated  by  very  sharp  resonant  pea)cs  whose  levels 
would  increase  without  bound  if  the  resolution  in  frequency  were 
increased.  This  is  because  of  the  absence  of  dissipation  in  the 
system.  In  Figure  3B  the  angular  separation  between  source  and 
receiver  is  the  maximum  of  180®.  The  feature  which  distinguishes 
this  plot  from  the  previous  one  is  the  increased  density  (factor 
of  2)  of  resonant  peaks,  simply  due  to  the  factor  cos(n0)  in  the 
summation  which  vanishes  for  odd  values  of  n  when  (;>=90®. 


5.  CONCLUDING  REMARKS 


The  Green's  functions  in  the  region  between  parallel  planes 
on  which  either  Neumann  or  Dirichlet  boundary  conditions  apply 
have  been  presented  in  a  form  which  facilitates: 

(i)  Theoretical  analysis  of  certain  boundary 
value  problems  in  finite  axisymmetric 
regions. 

(ii)  Numerical  evaluation  of  the  sound  fields 
without  incurring  problems  of  severe  ill- 
conditioning. 

A  closed-form  expression  for  the  sound  field  in  an  axisymmetric 
channel  and  some  numerical  results  have  been  obtained  as  an 
example.  The  availability  of  Green's  functions  result  in  their 
application  to  practical  problems.  For  example,  James  C9D  has 
used  the  Green's  functions  of  equation  (2.10)  to  investigate  the 
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